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Abstract 

Symbolic dynamics for homoclinic orbits in the two-dimensional symmetric map, 
Xn+i + cx n + = 3x^, is discussed. Above a critical c*, the system exhibits 

a fully-developed horse-shoe so that its global behavior is described by a complete 
ternary symbolic dynamics. The relative location of homoclinic orbits is determined 
by their sequences according to a simple rule, which can be used to numerically locate 
orbits in phase space. With the decrease of c, more and more pairs of homoclinic 
orbits collide and disappear. Forbidden zone in the symbolic space induced by the 
collision is discussed. 
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1 Introduction 



This paper concerns locating homoclinic orbits (HOs) for maps. The topic received some re- 
cent interests because it is closely related to finding spatially localized oscillatory solutions in 
one-dimensional lattices (see jl|] and references therein). Numerical methods usually involve a 
shooting strategy or an assumption of the initial topology of HOs (2|, ^, [|, |5|. In both cases, 
a deep understanding of the qualitative behavior of maps is important. As a coarse-grained 
description of dynamics, symbolic dynamics(SD) has been productively applied to one- and 
two-dimensional maps [p. The effectiveness of SD relies on the hyperbolic nature of the under- 
lying dynamics. Therefore, SD should be a powerful tool to describe the qualitative behavior 
of HOs when the system is dominated by unstable motion. In this paper, we shall discuss this 
through a specific map. 

The system we consider is the map with its orbits {x n } obeying the relation 

x n+1 + cx n + x n _i = 3x1, (1-1) 

which has been studied in H f7j. The coefficient 3 can be re-scaled to an arbitrary positive 
number, so only c is the parameter of dynamical significance. Let (a n ,b n ) = (x n ,x n -i), Eq.(l) 
can be rewritten as a two-dimensional conservative map, 

(a n+1 , b n+1 ) = F(a n , b n ) = (3a^ - ca n - b n , a n ). (1.2) 

The system contains a fixed point O = (0, 0), which is unstable when |c| > 2. As in Ref. @, we 
consider the orbits that are homoclinic to point O at c > 2. The paper is organized as follows. 
Section 2 gives a preliminary discussion of the map. Section 3 contains a general description of 
its SD. In Sec. 4, symmetric HOs are discussed in detail. This is followed by a brief summary. 



2 Preliminary 

2.1 Symmetries 

Map (1.2) is invariant under the spatial reversion (a, b) — > (—a, —b). Moreover, it has two inverse 
symmetries, I\ : (a, b) — > (b, a) and I2 : (a, b) — > (—6, —a). I\ and I2 satisfy I\FI\ = I2FI2 = F^ 1 
and l\ = l\ = I. The inverse symmetries determine four fundamental symmetry lines, 



2 



sh : Ii(a, b) = F(a, b) or 3a 3 - ca - 2b = 

sl 2 : h(a,b) = (a,b) a -b = , , 

sl 3 ■ h(a,b) =F(a,b) a = 1 ' ' 

sU : I 2 (a, b) = (a,b) a + b = 

For points on these lines, their forward and backward trajectories are related by simple rules. 

Specifically, F k (X) coincides with I 1 F 1 ~ k (X),I 1 F~ k (X),I 2 F 1 - k (X) or I 2 F' k (X) if X belongs 

to sli,sl 2 ,sls or s/4, respectively. Therefore, if X further belongs to the unstable manifold of a 

fixed point, it must be a homoclinic or herteroclinic point. This property was used to locate 

HOs with symmetry Q. 

2.2 Normal form for stable and unstable manifolds 

Numerical investigation of HOs relies essentially on a precise computation of the stable and 

unstable manifolds. For canonic systems, the invariant manifold of a unstable fixed point 

Xo = (ooi&o) can be determined by the Birkhoff normal form|J. The manifold is given by 

the parameterized curve X(t) = (a(t),b(t)) with 

00 00 
a(t) = ^ ai f and b{t)=^bit\ (2.2) 

i=0 i=0 

where coefficients {a,i,bi} are determined by requiring 

X(\t) = F(X(t)). (2.3) 

Here A is one of the eigenvalues of F at Xq, specifically, |A| > 1 for the unstable manifold and 
|A| < 1 for the stable one. Take Xq = O as an example. Its unstable manifold can be written as 

O u (t) = (G(Xt),G(t)), (2.4) 

where A = -(c + Vc 2 - 4)/2 and G(t) = git + g 2 t 2 + g-zt 2 "... with g\ = 1 and 

3 

9k = x k + X~ k + c ^ 9q9r9s (2 ' 5) 

q+r+s=k 

for k > 1. When k — > 00, gt drops very fast if c is not too small (e.g. c > 3). In our numerical 
study, O u {t) at t ~ 1 is computed by the series expansion with a reasonable cutoff. Applying 
F on a precisely determined piece of O u (t) (2.4), we may obtain O u (t) at large t with high 
precision. 
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3 Symbolic dynamics: general description 



3.1 Fully-developed horse-shoe and complete symbolic dynamics 

The global dynamics of (1.2) at large c is most conveniently accounted for based on the relation 
between the stable and unstable manifolds of the fixed points A = (xd, Xd) and B = (—Xd, —xa), 
where Xd = \/(c + 2)/3. Figure la shows these invariant manifolds at c = 5.7, denoted by 
A S ,A U ,B S and B u , respectively . The first part of A s , ACE, is approximately a horizontal line. 
F^ 1 maps AC to ACE and two additional nearly horizontal segments, EG and GJ, forming a 
s-turn. The first three parts of A u , B s and B u can be generated from A s by taking advantage 
of the symmetry. 

Under the action of F , the area confined by the first segments of A s , A u , B s and B u (the 
diamond-shaped ACBD in Fig. la) is transformed to the s-shaped region ACEGJBDFHK, 
indicating that the global dynamics exhibits a well-developed horse-shoe. This is schematically 
shown in Fig. lb. After taking into consideration of the symmetry, we can see that this geometric 
structure is guaranteed by the intersection between EG J consisting of the second and third 
segments of A s and AD, the first part of A u . The points of intersection correspond to two orbits 
that are homoclinic to point A. In fact, the two points lie on sl\ and hence correspond to HOs 
with reflection symmetry. With the decrease of c, the two points move toward each other and 
collide at c = c* ~ 5.453254835, leading a tangency of A s , A u and sl\. 

As a prototype of chaotic dynamics, the global behavior of horse-shoe is well-known. In our 
case, it is conveniently described by a complete ternary SD. 

Symbolic representation of points All stable (or unstable) manifolds, that generally refer 
to invariant vector field and shrink under the action of F (or i* 1-1 )^], can be naturally divided 
into three groups and coded with +, and — , respectively (see Fig. 1). Correspondingly, a 
point X in phase space is represented by a bilateral symbolic sequence 

S{X) = ...s-28-i • ao*i*2- = S-PO • S+(X), (3.1) 

where € {+,0,—} is the code for the stable (or unstable) manifold on which F k (X) (or 
F k+1 (X) ) is located. For example, the symbolic sequence for points C, J, D and K in Fig. la 
are -°° • +°°, -°° • -+°°, +°° • -°° and +°° • +-°°, respectively. 

The completeness of SD means that each sequence is realized by at least one phase space 
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point. By definition we have 

SoF(X) = FoS(X), (3.2) 

where T is the right shift of " • " with respect to the sequence. Therefore, a bilateral symbolic 
sequence without " • " can be used to represent an orbit without its initial point specified. For 
example, the symbolic sequences of the fixed points A, B and O are +°°,— 00 and O 00 , respectively. 

Order of symbolic sequences A point is indicated by its forward sequence 5+ (hence stable 
manifold) and backward sequence S- (unstable manifold). The location of stable manifolds in 
the phase space induce a nature order to forward symbolic sequences, •sosis^--- > •s' s' 1 s' 2 ... if 
the stable manifold of •SQS1S2... is above that of •s' s' 1 s' 2 .... The assignment of letters implies 

•+ > .0 > •-, (3.3) 

i.e., any sequence beginning with " + " is larger than any sequence beginning with "0", which 
is in turn larger than any sequence beginning with " — ". The comparison of two arbitrary 
sequences relies on the fact that + and — preserve the order while reverses it. (Note that the 
eigenvalues of fixed points A and B are positive while those of point O are negative.) Assume that 
S = siS2-.-s n s n+ i... and S' = siS2---s n s' n+1 ... and »s„+i > *s' n+1 , then »S > »S' (or »S < •S") 
if the number of "0" in siS2---S n is even (or odd). For example, the first three segments of A s 
correspond to •+°°, •0+°° and • — respectively. It can be easily verified that the first 
among the three is the largest forward sequence, and that there is no sequence which is smaller 
than the second and at the same time larger than the third. Similarly, backward sequences can 
be ordered in accordance with the location of their unstable manifolds in the phase space. We 
have 

...S3S2Si« > ...SgS^si* •S1S2S3... > •s' 1 s' 2 s' 3 .... (3.4) 

We shall indicate a point, an orbit or a piece of invariant manifold with the same symbolic 
sequence if there is no danger of causing confusion. 

3.2 Under-developed horse-shoe and truncated symbolic dynamics 

When c < c*, the second and third segments of A s do not cross the first segment of A u . In 
this case, there are some "forbidden" sequences, e.g. + oo 0+ oc and +°° — h°°, which cannot be 
realized by any real orbits. Consider a sequence S = S-S+. When c > c*, it corresponds to a 
real orbit on which there is a point being the intersection of stable manifold •S + and unstable 
manifold 5_«. When we reduce c, the two segmens of manifold will vary continuously with c 
and they must be tangent at a critical value of c = q, before S is forbidden. The tangency can 
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be viewed as a collision of two orbits. Since the two colliding orbits are infinitely close when 
c — > Cfe, generally their corresponding sequences differ by only one symbol, say, one is S1OS2 and 
the other is Si — S%. (A logically more meaningful statement is that, to maintain the simple 
order of invariant manifolds, the partition line that separates SiO* and Si — • should pass the 
tangent point ||]). The tangency of stable manifold «S2 with unstable manifolds SiO» and Si — • 
implies a forbidden zone in the space of sequences or symbolic space||]. Specifically speaking, 
sequence S' = S[ • S' 2 with 

Si - • < S'i» < SiO • and • S' 2 > »S 2 (3.5a) 

corresponds to no point in the phase space (Fig. 2a). Similarly, if the two colliding orbits are 
S1OS2 and Si + S2, the forbidden zone is given by 

Si + • > Si* > SiO • and • S' 2 < *S 2 . (3.56) 

Of course, any sequence S' with certain T k {S') satisfies (3.5) is also forbidden. 

Due to symmetry, it is not exceptional in our case that the sequences of the two colliding 
orbits differ by two symbols. The forbidden zone due to such collision will be also discussed. 

4 Symmetric homoclinic orbits 

We now focus on the orbits that are homoclinic to the fixed point O. The unstable manifold 
of O is given by (2.4). From the symmetry, the stable manifold can be written as O s (t) = 
(G(t),G(\t)). By definition, HOs correspond to the solutions of equations G(Xt) = G(t') and 
G(t) = G(Xt'). However, the equations have little practical use because G(t) changes violently at 
large t. The more useful knowledge is the global location of the stable and unstable manifolds. 
When c > c*, each sequence O^SO 00 corresponds to a HO. To locate the orbit in the phase 
space, we can scan initial point X on 0°°», the first segment of O u , searching for the target 
S_|_(X) = SO 00 . For a symmetric HO, the target can be replaced by F k {X) £ sli or alternatively 
the scan can be performed on a symmetry line. In all cases, the order of forward sequences can 
be used for an effective searching strategy. Having a HO at large c, we can trace it to small c 
until it collides with another HO and disappears. We shall not go into details about the locating 
of HOs. We shall focus on the compatibility between HOs at c < c*. 

In the following discussion, the stable and unstable manifolds are always referred to as O s 
and O u . We first give another description of the forbidden zone induced from the tangency of 
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stable and unstable manifolds, or the collision of two HOs. Let X and X' be two homoclinic 
points which collide each other when c — ► q>. For c > c*, the stable and unstable manifolds 
connecting X and X' form a closed curve C?(X, X') (see Fig. 2), which shrinks to a single point 
when c = c&. Note that both O s and O u are continuous and self-avoiding curves. Consider a 
homoclinic point enclosed by 0(X,X r ) or lying on X'). Being the intersection of stable 

and unstable manifolds, the point cannot cross 0(X,X') when c is continuously varied. The 
area enclosed by 0(X,X') at c > c* vanish before c = c&. Therefore, C(X, X') defines the 
forbidden zone corresponds to the collision of X and X' . 

It can be seen that the forbidden zone defined by 0(X, X') coincides with that defined by (3.5) 
if S(X) and S(X') differ by only one symbol (see Fig. 2b). To study the evolution of symmetric 
HOs, we need also consider the the case that S(X) and S(X') differ by two symbols. For this 
purpose, the structure of stable and unstable manifolds should be studied in detail. Consider 
the 3 fc segments of stable manifold of the form •SO 00 , where S exhausts all the ternary strings 
with length \S\ = k. The relative locations of these segments of manifold in the phase space 
are determined by the order of »S. We ask that how the 3 fc segments are connected to form a 
continuous curve. All these segments are generated from "O 00 by F~ k , i.e., J 7 ~ k (S»0 oc ) = •50°°. 
We may cut the segment "O 00 , according to backward sequences, into 3 k pieces, each of which 
are coded with SV The way to joint these 3 k pieces into the whole segment "O 00 is determined 
by the order of the 3 fc strings S». After expanding under F~ k , the segment "O 00 becomes the 
3 fc segments of "SO 00 , so they should be connected also according to SV In the simplest case 
of k = 1, three segments of O s from top to bottom are • + 0°°, "O 00 and • — 0°°. Forming an 
s-turn, • + 0°° and • — 0°° joint with "O 00 at its left and right ends, respectively. For k = 2, the 
nine segments of O s are ordered from top to bottom as 



as shown in Fig. 3 schematically by horizontal lines. In the figure S is also taken as the abbre- 
viation for •S0°°. These segments are connected end to end according to the following order of 
backward SV 



The arrangement is shown in Fig.3. We see that pair (00, +0) joint at their left ends, (+0, H — ) 
joint at their right ends, (H — , 0— ) joint at their left ends, and so on. In general, the sequence 
pair of connected segments of stable manifold take one of the following forms, 



++, +0, +-, 0-, 00, 0+, -+, -0 and -- 



++, 0+, -+, -0, 00, +0, +-, 0- and 




(i) S = - m 0S and S' = - 
ii) S = + m 0S and S' = + : 



+ S , 
- So, 



7 



where m > 0. Segment pair of type (i) are connected at their left ends while pair of type (ii) 
are connected at their right ends. By taking advantage of the symmetry, a parallel description 
can be given to the unstable manifold. 

Now we study the forbidden zone implied by the collision of two symmetric HOs. Let us first 
consider two examples. In the first case, S(X) = 0°° • + + +0°° and S(X') = 0°° - • + + + -0°°. 
From the above discussion, 0°°» and 0°° — • are connected at their upper ends while • + + + 0°° 
and • + + H — 0°° are connected at their right ends. The forbidden zone enclosed by 0(X, X') is 
shown in Fig. 2c. We see that it is the union of two forbidden zones, one is defined by 0(X, X") 
and the other by 0(X' , X"), where S(X") = 0°° • + + + - 0°° (or 0°° - • + + + 0°°). 

The next example is S(X) = 0°° - • 00 - 0°° and S(X) = 0°° - + • + Now 0°° - 0* 
and 0°° — h* are connected at the bottom. On the other hand, »00 — 0°° and »0 H — 0°° are not 
joint neighboring segments on O s . They are connected by • — — 0°° and • — I — 0°°. The area 
enclosed by 0(X, X') is shown in Fig. 2d, which again turns out to be the union of two regions 
enclosed by 0(X, X") and 0(X', X") with S(X") = 0°° - • + -0°° (or 0°° - + • 00 - 0°°). 

In general, suppose S{X) = 0°°5i0 • S 2 QS^°° and S(X') = 0°°5iai • S 2 a 2 S 3 Q°°, where 
^1,2 G {+,—}. This covers all the colliding symmetric HOs with two different letters. As 
seen in the two examples, it can be shown that the forbidden zone enclosed by 0(X, X') is 
the union of that defined by 0(X,X") and 0{X',X"), where S(X") = 0°°5i0 • S 2 a 2 S 3 0°° (or 
0°°Siai • S 2 0S 3 0°°). In other words, when forming the forbidden zone, the collision of X and 
X' can be regarded as the triple collision involving X, X' and X" . Since either the pair of S(X) 
and S(X") or the pair of S(X') and S(X") differ by one symbol, the corresponding forbidden 
zone is given by (3.5). Thus, the rule can be used to construct the forbidden zone for the above 
general case. 

We have numerically calculated all symmetric HOs O 00 ^ 00 with \ S\ < 7 and traced them from 
c = 6 to collision parameters q>. Among the 160 orbits, only 0°° ± 0°° and 0°° ± +0°° persist till 
c = 2. We obtained 94 collisions, each of them corresponds to a tangency of O s and a symmetry 
line (see Fig. 4 for some examples). Due to symmetry, only half of the collisions are independent. 
They are ordered with descending q, and listed in Tab. 1. We further examine the compatibility 
between symmetric HOs. Collision of a pair of symmetric HOs may implies the elimination of 
some others. For example, it can be verified that the forbidden zone of the colliding HO pair 

(S, S') = (H 1 — , +00—) includes the HO sequence pair ( — I — , —0—). Therefore, the collision 

of HO pair ( — I — , —0—) should occurs at a larger c than that for (H 1 — , +00—). In fact, the 

former appears as the 9-th collision in Tab. 1 while the latter is the 16-th. We find that all such 
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restrictions in the 47 pairs of sequences in Tab.l are in accordance with the order of q, (Fig.5). 



5 Summary 

In the above we have discussed the symbolic dynamics of system (1.2) with a focus on the orbits 
that are homoclinic to point O = (0,0). When c > c*, the global dynamics is described by a 
complete ternary symbolic dynamics, so that each sequence O^SO 00 corresponds a homoclinic 
orbit. The location of these orbits are qualitatively determined by their corresponding sequences 
according to simple rules. With the decrease of c, more and more symmetric homoclinic orbits 
collide in pair and disappear. Based on the geometrical constraint that any homoclinic point 
cannot cross a region enclosed by the stable and unstable manifolds, the compatibility between 
homoclinic orbits is discussed. 
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Figure Captions 



Fig.l (a) First three segments of A s , A u , B s and B u at c = 5.7. The dotted lines (a, b = ±y/c/3) 
give a convenient partition of manifolds. (b)Schematic diagram of hose-shoe. The square ACBD 
is mapped by F" 1 to the shaded area. 

Fig.2 Schematic diagrams to show forbidden sequences due to tangency between the stable and 
unstable manifolds. Figure (a) explains rule (3.5a). The shadowed areas in (b-d) represent the 
forbidden zones when X and X' collide. 

Fig. 3 Arrangement of the first 9 segments of O s . 

Fig.4 Tipical tangencys of O u (thick line) and symmetry line (thin line). Each tangeny represents 
a collision of symmetric homoclinic orbit pair. 

Fig.5 Compatibility between symmetric homoclinic orbit pairs listed in Tab.l. A dot at means 
that the collision of i-th pair implies the elimination the j-th pair. Note that all dots are located 
within the region i > j. 



10 



Tab.l Critical parameter of symmetric HOs 
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